Appl. No. 10/062,483 

REMARKS 

Status of the Claims 

Claims 1-20 are pending in this application. Claims 15-20 have 
been withdrawn from consideration. No claims have been canceled, 
added or amended. 

Rejection under 35 USC 112, first paragraph 

The Examiner rejects claims 1-14 as not enabled by the 
specification. The Examiner states that all compounds that bind to 
the SSI and SS2 subunits of a sodium channel are not enabled by the 
specification. The Examiner states that only tetrodotoxin, 
saxitoxin and derivatives thereof are enabled. Applicants traverse 
the rejection and respectfully request the withdrawal thereof. 

Applicants submit that the Examiner has failed to meet the 
burden of presenting a prima facie case as to why the claims would 
not be enabled. See In re Wright, 27 USPQ2d 1510 (Fed. Cir. 1993) . 
Wright, citing In re Marzocchi, 169 USPQ 367, 369 (CCPA 1971) 
states 

When rejecting a claim under the enablement requirement of 
section 112, the PTO bears an initial burden of setting forth 
a reasonable explanation as to why it believes that the scope 
of protection provided by that claim is not adequately enabled 
by the description of the invention provided in the 
specification of the application; this includes, of course, 
providing sufficient reasons for doubting any assertions in 
the specification as to the scope of enablement. If the PTO 
meets this burden, the burden then shifts to the applicant to 
provide suitable proofs indicating that the specification is 
indeed enabling. 
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The Examiner has failed to meet this initial burden. 

Even if the Examiner had met this burden, Applicants have 
submit that the present specification provides adequate 
enablement. It is not necessary that the specification provide 
examples of every possible operable embodiment, only that the 
specification, taken with what was known at the time the 
application was filed, describe to the skilled artisan how 
embodiments throughout the scope of the claims can be made 
without undue experimentation. 

The present specification describes a number of compounds 
effective in the present invention, and this is acknowledged by 
the Examiner. Furthermore, Applicants also submit herewith 
evidence that establishes that the state of the art at the time 
the application was filed was such that certain structural and 
functional similarities between compounds that bind SSI and SS2 
subunits were recognized. Applicants submit that one of 
ordinary skill in the art would be able to identify compounds 
that bind SSI and SS2 subunits for use in the present invention. 
Absent a preponderance of evidence from the Examiner that the 
skilled artisan could not identify compounds that bind to SSI 
and SS2 subunits that work in the present invention, a 
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presumption arises that the full scope of the claimed invention 
is enabled by the specification. 

Applicants submit herewith three journal articles: Yang, L. 
et al. "Actions of Chiriquitoxin on Frog Skeletal Muscle Fibers and 
Implications for the Tetrodotoxin/Saxitoxin Receptor" , Journal of 
General Physiology , Vol. 100, pages 609 to 622, published October 
1992; Fozzard, H. et al . , "The Guanidinium Toxin Binding Site on 
the Sodium Channel", Japanese Heart Journal , Vol. 37, No. 5, pages 
683 to 692, published September 1996; and Lipkind, G. et al . , "A 
Structural Model of the Tetrodotoxin and Saxitoxin Binding Site of 
the Na+ Channel", Biophysical Journal , Vol. 66., pages 1 to 13, 
published January 1994. Each of these references supports 
Applicants' position that the claims are enabled. 

For example Lipkind et al . describes a three dimensional model 
of the SSI and SS2 subunits with TTX bound thereto. The model is 
so detailed as to show contacts of the toxin with the protein. 
Such a detailed model allows the skilled artisan to predict the 
structure of an effective sodium channel blocker. The authors 
state in the introduction that TTX and STX have virtually identical 
biological activity due to a common guanidinium group that is 
positively charged at physiological pH. The authors note that a 
guanidinium group is common to compounds that bind SSI and SS2 
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subunits of the sodium channel and is important for blocking the 
channel . 

See also the Abstract of Yang et al . , which discloses that 
chiquirotoxin (CqTX) , another sodium channel blocker, is only 
different from TTX in the substitution of a glycine residue in 
place of a methylene hydrogen of the C-ll hydroxyl methyl function. 

Lastly, Fozzard et al . explains that guanidinium is a common 
feature of compounds that bind to the sodium channel . 

Applicants submit that compounds that bind to SSI and SS2 
subunits of the sodium channel are easily determined without undue 
experimentation, particularly given knowledge of a chemical 
structure that is a feature of such compounds. A guanidinium group 
is likely to be present in an effective compound. This knowledge 
provides a starting point for any expected experimentation. To 
this knowledge, Yang et al . add a technique for determining 
effectiveness of the compounds selected. Yang et al. use the patch 
clamp method to measure the channel current in the presence and 
absence of an inhibitor. 

The state of the art at the time the application was filed was 
thus such that the skilled artisan could predict the structure of a 
SSI or SS2 -binding sodium channel blocker and could test that 
compound for effectiveness using assays known in the art. Thus, no 
undue experimentation is required to practice the present invention 
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throughout the scope of the present claims. Therefore, Applicants 
submit that the claims are enabled and the instant rejection should 
be withdrawn. 

Conclusion 

As Applicants have addressed and overcome the rejection in 
the Office Action, Applicants respectfully request that the 
rejection be withdrawn and that the claims be allowed. 

Should there be any outstanding matters that need to be 
resolved in the present application, the Examiner is respectfully 
requested to contact Kecia Reynolds (Reg. No. 47,021) at the 
telephone number of the undersigned below, to conduct an interview 
in an effort to expedite prosecution in connection with the present 
application. 
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If necessary, the Commissioner is hereby authorized in this, 
concurrent, and future replies, to charge payment or credit any 
overpayment to Deposit Account No. 02-2448 for any additional fees 
required under 37 C.F.R. §§ 1.16 or 1.17; particularly, extension of 
time fees. 



Respectfully submitted, 



BIRCH, STEWART, KOLASCH & BIRCH, LLP 



By. 



Mark 



JVNuell, #36,623 



DRN/KJR/jao 



P.O. Box 747 

Falls Church, VA 22040-0747 
(703) 205-8000 
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ABSTRACT Biophysical evidence has placed the binding site for the naturally occurring marine toxins tetrodotoxin (TTX) 
and saxitoxin (STX) In the external mouth of the Na" channel ion permeation pathway. We developed a molecular model of 
the binding pocket for TTX and STX, composed of antiparallel 0-hairplns formed from peptide segments of the four S5-S6 
loops of the voltage-gated Na + channel. For TTX the guanidinium moiety formed 3alt bridges with three carboxyls, while two 
toxin hydroxyts (C9-OH and ClO-OH) interacted with a fourth carboxyi on repeats I and II. This alignment also resulted in a 
hydrophobic interaction with an aromatic ring of phenylalanine or tyrosine residues for the brainli and skeletal Na* channel 
isoforms. but not with the cysteine found in the cardiac isoform. In comparison to TTX, there was an additional interaction 
site for STX through its second guanidinium group with a carboxyi on repeat IV. This model satisfactorily reproduced the ef- 
fects of mutations in the S5-S6 regions and the differences in affinity by various toxin analogs. However, this model differed 
in important ways from previously published models for the outer vestibule and the selectivity region of the Na"** channel 
pore. Removal of the toxins from the pocket formed by the four 0-hairpins revealed a structure resembling a funnel that ter- 
minated in a narrowed region suitable as a candidate for the selectivity filter of the channel. This region contained two car- 
boxyls (Asp 304 and Glu** 2 ) that substituted for molecules of water from the hydrated Na* ion. Simulation of mutations in this 
region that have produced Ca 2 ^ permeation of the Na+ channel created a site with three carboxyls (Asp 38 *, GIu** 2 , and 
Glu ,71d ) in proximity. 



INTRODUCTION 

Tetrodptoxin (TTX) and saxitoxin (STX) arc naturally oc- 
curri ng high affin ity marinelQ xi^ 
volta^sfinsitia^e- Na^ change! rHjlle7l992TCaneraU. 1992)' 
Their overlapping binding sites are most probably located in 
the outside mouth of the ion permeation pathway. Although 
the cjifi mjca] stn icjures oJTTX and STX (see Fig. 1) aregu ite 
4\\ffRTrj\t i rt^'y hjHpfli cal acti vi ties ar c virtually idejiticaL 
partly as a consequence of a share d guanipligjum group that 
is pj2sjthgly ctyged atph ysioIg^ica LpH. One or more nega- 
tively ch^ge^^arboxyigroupS on the channel protein have 
been shown to play an essential role in toxin binding by either 
pH titration or by treatment with caiboxyl-modifying re- 
agents (Hille, 1968; Shrager and Profera, 1973). The gua^- 
^n jdiniu m Jon itself is jj g gneani to the Na* channe l (Hille, 
1971), and bothpcmieant and non-penneant (blocking) cat- 
ions can compete for the toxin site (Narahashi, 1974). On the 
basis of such data, Hille (1975) suggested that the gua- 
nidinium group of TTX is bound to a rectangular ring of six 
oxygens located in the walls of the channel pore at its se- 
lectivity filter, with the cationic guanidinium forming a salt 
bridge to one negatively charged acid group and with other 
polar groups of the toxin interacting with the remaining oxy- 
gens so as to occlude the pore. 
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In I2£4J^uma and coworkers (Noda et al. f 1984) cloned 
the first of a fa mily of isofor ms of the principal Sub uniTo? 
the Na* channel These Na*~ channel Ta-subunits^reTarge 
glycoproteins of about 2000 amino acids, organized into four 
highly homologous domains or repeats (I-IV). Each domain 
includes six transmembrane segments (S1-S6) with high 
a-helical potential (Numa and Noda, 1986). On the basis of 
biophysical and chemical data the extracellular segment be- 
tween S5 and S6 of each repeal is thought to fold back into 
the membrane to form part of the ion-conducting pathway 
(Guy and Conti, 1990). A series of point mutations in this 
S5-S6 segment by Stuhmer, Numa, and colleagues and oth- 
ers (Nodaet aL, 1989; Terlau et aL, 1991; Satin et aL, 1992; 
Backx et aL, 1992; Heinemann et aL, 1992a; Kontis and 
Goldin, 1993) has identified several specific amino acids 
essential for toxin block and for isoform differences in toxin 
affinity. 

The size and solubility characteristics of the Na* channel 
preclude determination of the three-dimensional structure of 
mis large molecule at this time. However, the primary struc- 
tures of several isoforms are known, along with effects of 
point mutations that affect toxin binding. We used this in- 
formation to develop a model of the toxin binding pocket 
from parts of the S5-S6 loops that have been identified as 
probable parts of the outer vestibule of the channel pore. The 
model was based on the simple question if a reasonable struc- 
tural organization of the amino acids in this region could be 
arranged around the toxin molecules so as to form an. inter- 
action that satisfied the available analog and mutation data. 
Working with only a small fraction of the entire protein se- 
quence, we found a predicted structure of these short peptide 
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segments that interacted well with the toxins, yielding a bind- 
ing pocket conforming to the available experimental toxin 
analog and channel point mutation data and to channel iso- 
fonn differences in toxin block. 



EXPERIMENTAL RESULTS USED IN MODEL 
DEVELOPMENT 

The construction and constraints of this model depend 
strongly on experimental studies of block by toxins and their 
analogs, and on the available data on point mutations in 
cloned Na"*" channels. 

Blockade by various guanidinium toxins 

Kao and coworkers (see review by Kao, 1986) have iden- 
tified the active groups in TTX and STX. In the case of TTX 
these are the guanidinium group and the hydroxy Is at C9 and 
CIO (see Fig. 1). The guanidinium group is thought to form 
an ion-pair with an anionic site of the receptor, while the 
(g^QH ^ndClO-OH formhydrog en bonds with other sites 
(Kao anoWalker, 1982)nrhe latter suggestion waTcon- 
^ firmed by the low (0.01-0.02 X) affinity of the analog an- 
hydroretrodotoxin, where C4 and C9 are joined by an oxygen 
bridge (NarahasJbi et aL, 1967; Kao and Yasumoto, 1985). 
4 ? Epitetrojo toxin < in which the positions of the -H and -OH 
were reversed from TTX, 4s about 0.5 X as p otent as TTX 
Therefore, the change in activity of anhydrotefroHotoxin can 
be explained as loss of the C9-OH group and one of the 
hydrogen bonds to the binding site. Modifications of both 
C9-OH and ClO-OH in the case of tetrodonic acid led to a 
totally ineffective molecule (Narahasbi et al, 1967). 
6Jjipi tetrodotoxin and 11-deoxvtetrodotoxin h ave low bjnj- 
in^ju!finiry.(about 0.01 x) to^Na^d^nn^ (Yang et al., 
1992a). However, nortetrodotoxin and nortetrodotoxin al- 
cohol, which cannot form corresponding hydrogen bonds, 
were somewhat less potent than TTX (Kao, 1982), implying 
that the C6 end of TTX is not very important for channel 
block. 

In prior studies of STX, the T^&^gianidirimm flmup and 
♦ two C13 , hydiQsyjs. wercjgejarjged^ nl Joj^SXX 
'^cjiig (Scrichanz, 1984; Kao, 1986). The two hydroxyls 
can be reduced selectively to produce cither the a or the p 
epimer of saxitoxinoL which have blocking affinity of only 
0.01 X (Koehn et aL, 1981; Kao et aL, 1985). Koehn et aL 
(1981) concluded thai the carbamoyl group of STX (H2NCO- 
group as substitute at C6) contributes to, but is not essential 
for channel blockade, because decarbamoylsaxitoxin (dc- 
STX) is about 0.2 X as active as STX Besides, an acetyl 
derivative of deSTX, in which an amino group was substi- 
tuted by a methyl group, has shown the same activity (Mahar 
et al., 1991). Nevertheless, the potency of decarbamoylox- 
ysaxitoxin is only O.0O8 X that of STX (Yang et al., 1992b). 

Yang and Kao (1992) identify seven possible interaction 
sites in the structures of TTX and STX that they conclude are 
important, but a more conservative analysis suggests that 
only the^guanidtnlum group and the two hydroxyls a t C9and 
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£XO Jfot TTX) or C12 (f r STX) are the most criticalsites. 
These sites show remarkable stereospecific similarities and 
appear to determine the high toxin specificity. When the gua- 
nidinium group of TTX and the 7,8,9-guanidiniuro group of 
STX are aligned, then the hydroxyls at C9 and C10 for TTX 
and the rwo hydroxyls at C12 for STX are also closely 
aligned (Kao and Walker, 1982). The positive charge in the 
second 1 ,2,3-guanidinium group of STX may also be in- 
volved, but direct evidence is lacking. Finally, chiriquitoxin 
(CqTX) differs from TTX only by a glycine residue in place 
of a hydrogen of the hydroxymethyl group at position Cll, 
but at pH 7.25 it has equal block potency as TTX (Yang and 
Kao, 1992). Since glycine has high potency for formation of 
hydrogen bonds but its presence does not modify binding, the 
CI 1 site is most probably outside of the vestibule and not able 
to interact with the binding pocket. 

Inference for binding from mutation studies 

Site-directed mutations of , the conserved amino acids of the 
rat brain n (rBrll), the adult skeletal, and the cardiac 
or-subunits of the Na* channel have provided important in- 
formation about the amino acids residues required for TTX 
and STX binding. We are indebted to the extensive work of 
Stuhmer, Numa, and colleagues for explorations of the amino 
acids previously suggested by Guy and Conti (1990) to be 
part of the pore and named SS1-SS2 or the P loop. A sing le. 
point mutation j)f glutamic acid 387 (rBrll num ber) to glu- 
lamine (S 87Q) results Jn^c omplete loss of T^andjTX 
jtlock (Noda^ 5^989). Glu 3 ^ belongs To the SS2 segment 
of repeat I (FigT^r Terlau et al. (1991) made systematic 
mutations in residues of the SS2 segments of repeats MV. 
Neutralization of any of the six conserved residues of Glu or 
Asp reduced the sensitivity to TTX and STX block by at least 
three orders of magnitude (Table 1). Charge mutations at 
other adjacent positions produced smaller changes in toxin 
block. Terlau et al. (1991) suggested that there are two clus- 
ters of negatively charged residues at equivalent positions in 
the SS2 segments of repeats I-IV (marked by arrows in Fig. 
2). The most dramatic effect on toxin block in their study was 
produced by neutralization of the carboxyl residues on re- 
peats I and II. Neutralization of aspartic acid at position 1717 
of repeal TV (D1717N) abolished STX block, while partially 



TABLE 1 Effects of single mutations on TTX and STX 
sensitivities and afnglo channel conductance 



Mutant 


Repeal 


IC50.TTX 


IQjo, STX 


Conductance 






nM 


nM 


pS 


Wad-type 




18 2 4 


1.2 ± 0.2 


15.4 


D384N 


I 


>ioooo 


>1000 


<0.1 


E3S7Q 




> 10000 


>1000 


3.1 


E542Q 


n 


>10000 


>1000 


<OJ 


E945Q 




2300 ± 160 


>1000 


8^ . 


D1426N 


m 


30 ±10 


8.9 ± 0.4 


14.0 


D1717N 


IV 


350 2 90 


>1000 


8.8 



Dau from Terlau tt al (1991). 
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reducing the block by TTX. In repeat III neutralization of 
Asp 1426 (mutant D1426N) had little effect on either toxin. 
Finally, the cardiac isofonn has a lower affinity for both 
toxins, which appears to be a consequence of the absence of 
au aromatic ring residue (Tyr or Phe) at position 385 in the 
rBril number system (Satin et al. t 1992; Backx et al., 1992; 
Keinemann et aJU 1992a). 

STRUCTURAL MOTIF FOR THE TTX/STX 
BINDING SITE 

For this model we infer that a residue is involved in a non- 
bonding or electrostatic interaction with the toxin if its sub- 
stitution changes the toxin blocking affinity 100-1000 X. 
Effects on binding affinity less thai 10 X can be physiologi- 
cally interesting, but beyond the accuracy of the present 
model. 

Scheme of principal sites of interaction 

In agreement with the experimental results of Tertau et al. 
(1991), we assumed that similar structural elements of the 
two toxins interact with residues of repeats I and II (Table 
1). The guanidinium group and the hydroxy Is C9-OH and 
C10-OH of TTX and the 7,8,9-guanidinium group and the 
bydroxyls at C12 of STX interact primarily with carboxyls 
at positions 384, 387, 942, and 945. The 1,2,3-guanidmium 
group of STX may interact with the carboxyl at position 1717 
of repeat IV, which has much less influence on TTX block 
(Tabie I). We propose that three of the carboxyls of repeats 
I and TI interact with the guanidinium group, while the fourth 
negatively charged residue forms two hydrogen bonds with 
neighboring hydroxy 1 groups of C9 and CIO for TTX and 
C12 for STX. Among the four carboxyl sites Glu 945 appears 
to be less important for TTX binding than the others, since 
the mutation E945Q reduces TTX sensitivity less than the 
other mutations. Most probably, Glu* 4,5 interacts with the 
TTX hydroxyls. Indeed, the hydrogen bond formed by the 
CI0-OH group is rather weak, as indicated by its easy de- 
protonizarion (Kao, 1986). On the other hand, both hydrogen 
bonds formed by the gem-diol group of STX are strong, so 
that Glu 945 is as important as the other carboxyls for STX 
binding (Konris and Goldin, 1993). It seems plausible that 
Glu 945 interacts with the toxin hydroxyls, and the other three 
carboxyls at positions 334, 387, and 942 form ion-pair salt 
bridges with the guanidinium group. 

We will consider an arrangement for toxin binding that 
includes four loci for TTX and five loci for STX (Fig. 1 ). This 
scheme is in accordance with the three resonance forms of 
the guanidinium, in which the positive charge is distributed 
between the N — H groups. Moreover, there are examples in 
proteins where different N — H groups of Arg form salt 
bridges with some carboxyl groups of Asp or Glu (Proegman 
et aL. 1978). Consequently, each of the three carboxyl resi- 
dues would contribute almost equally to interaction with the 
guanidinium group, and the fourth carboxyl group gives ap- 
proximately similar energetic contribution. Fersht et al. 



P. 3 

3 



Repeat II 

111 Mf 




QlFWgN 
Repeat I 



Repeat II Repeat IV 

MS MB X7X7 




Repeat I 

FIGURE 1 Molecular structures of TTX and STX and suggested amino 
acid interactions of the Na * rhnnnfl (see Fig. 2). There are four loci for TTX 
and five Jod for SIX. The toxin Structures flie from HOle (1992). 

(1985) showed that the energy of a salt bridge is approxi- 
mately equal to that of two neutral hydrogen bonds. 

Conformation of the SS1 and SS2 peptide 
fragments 

The first step in predicting a satisfactory protein structure for 
toxin interaction was to examine the probable structural 
properties of the peptide chains comprising the SSI and SS2 
segments from the four ( repeats of rBril (Fig. 2), talcing into 

v w 

F W 1 N L Y 

wiItmw 

Repeat III LtiQVATFKGWMglM 

1717 

Repeat IV LFQITTSAGWgGLL 
■ SSI 1 i SS2 1 

FIGURE 2 The SS1-SS2 regions (P-loops) qf repeats I-W of Lfle rat brain 
rj Na* channel. Positions of -certain critical mutations are marked (from 
Terlau et aL (1991)). 
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account the probably extracellular location of these seg- 
ments. SSI and SS2 contain many residues that have high 
probability for appearance in 0-strands or 0-structures: Val, 
Ue, Tyr, Cys,Trp, Phe, Thr, Met, and Leu (Creighton, 1993). 
For example, in the case of the SSI fragment of repeat II there 
are five bulky hydrophobic residues in the C-tennina] six- 
member segment; lie, Val, Leu, and Phe (Fig. 2). On the oiher 
hand, residues at the border between SSI and SS2 have a 
conformational preference to participate in reverse 0-tums: 
Gly, Asp, Ser, and Lys (Wilraot and Thornton, 1988). There- 
fore, this region has a high probability to form 0-hairpins, in 
which adjacent p-strarids are linked through /3-turns and pro- 
duce a structure that might exist at an interface between the 
hydrophilic surface and me hydrophobic core or the protein. 

In order to select liJcely positions for the /3-ruras in the 
peptide chains we used the procedure proposed by Chou and 
Fasman (197$), but taking into account newer values for the 
conformational preferences (P t ) fonnd on the basis of x-ray 
data of proteins (Williams et aL, 19S7). The values of A 
characterize the relative tendencies of residues to be involved 
in 0-turns (Table 2). Values of P x > 1.0 characterize amino 
acid residues favoring p-turns. According to Chou and Fas- 
man (1978), a j3-turn is predicted if a four-residue sequence 
has an average value </>) > 1 (and <i>) > (p^). Further, 
values of P ( for the two central residues of the j3-turn should 
be greater than those for the extreme residues. These con- 
ditions are reasonably met on the central portions of the four 
peptide fragments underlined in Table 2 ((PJ » 0,96, 0.99, 
1.02, and 1.12, respectively). The sequences preceding the 
proposed £-turns are composed of residues with low prob- 
ability of participating in turns. For repeat IV a second pos- 
sible rum region is seen, but we selected the first one for 
reasons presented later. In this manner, we chose the four 
^-hairpins of repeats I-TV, composed of two /3-strands con- 
nected in the center by 0-tums (as marked in Table 2), and 
consisting of 10 residues each. It was important to include 
all of the SS2 amino acids, in order to stabilize the 0-haiipin. 
The 1-2 adjacent residues on both ends also have a high 
probability of being in j3-strands, and the usual length of such 
/3-hairpins is 8-16 amino acids (Sibanda et al., 1989). Since 
these are no longer identical with the original definition of 



the SSI and SS2 segments, we will call them betal and beta2 
on the N and C ends, respectively. 

The conformations of the /3-hairpins are important for our 
modelled structure, because they determine the characteristic 
orientation of side chains (Fig. 3). We adopted the usual 
numbering for the four residues that comprise /3-tums (r, I + 
1, i + 2, and i + 3) (Creighton, 1993). Neighboring residues 
of the j3-strands alternate relative to the plane of the /3-hairpin 
(up-down-up...). Residues i and i 4- 3 belong to both the 
£-rura and the 0-strand. If residues i + 2, i + 3, and i -+■ 5 
of beta2 are located on the outside plane of the 0-hairpin, the 
residues i + 1, / + 4, and i + 6 are found on the opposite 
side. In repeat I we know from mutation studies that Asp 384 
(/ + 2) is important for toxin block. Therefore, Phe 385 (i -+- 
3), and Glu 387 (£ -h 5) would be important for toxin block, 
also consistent with mutation data, and other residues of the 
fragment that face away from the toxin site would be rela- 
tively unimportant In like manner we can predict for repeat 
II that only residues in positions j + 2> £ + 3, and i + 5 
(Glu 943 , Trp* 43 , and Glu 945 ) have the correct location to in- 
teract with toxin. Effects of these residues would be large 
(two to four orders of magnitude), while other residues would 
produce only small effects on toxin interaction by secondary 
changes in the binding structure. 

TTX interacted with rmmarily the segments of repeats I 
and II. These two repeats can be aligned to form a structure 
consisting of two antiparallel ^-hairpins adjacent to each 
other. This structure produces a cavity for binding of TTX, 
with the inner part consisting of six side chains of residues 
Asp 384 , Phe 385 , Glu 387 , Glu 942 , Tip 943 , and Glu 943 (Fig. 4). 
Presumably the /3-hairpins of repeats HI and IV form mainly 
the outer walls of the binding cavity for TTX, with the 
0-hairpins also arranged in an antiparallel fashion. For repeat 
m this results in the possibility of toxin interaction with the 
residues of Lys 1422 (i + 1), Met 1423 (/ + 4), and lie 1427 (i + 
6). The residues of interest for repeat IV in the arrangement 
are Ala 1714 (i + 1), Asp 1717 (/ + 4), and Leu 1719 (i + 6). 
Indeed, alignment of the ^hairpin of repeat IV in an anti- 
parallel relationship to that of repeat I on the front wall 
changes the directions of its side chains in space, -so that the 
side chains of residues / + 1, / + 4, and i + 6 are inside the 



TABLE 5 Conformational preferences Pf for residues In SS1^SS2 segments of the muscle channel and positions of B-tums 
(underlined) ^ 
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0.59 


R 
0,90 


V 
0.41 


L 
0.S7 


C 


G 
L2Z 


942 
E 
LQ1 


W 


I 

0.47 


945 
£ 
1.01 


T 

0.90 


M 
052 


Repeat LO 
P, 


L 

0-S7 


L 

U57 


O 
0.84 


V 
0.41 


A 


T 
0.90 


F 
0.59 


1422 
K 

LQ7 


G 

\m 


W 

0.65 


1425 
M 
0.52 


D 

1.24 


I 

0.47 


Repeat TV 


L 

0,57 


F 
0.59 


Q 
0.84 


I 

0.47 


T 

0.90 


T 
0.90 


s 

1-22 


1714 
A 

0.82 


G 
U7_ 


W 

0,65 


1717 
D 
L24 


G 

1.77 


L 
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• From Williams et &L (1987). 
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FIGURE 4 Structural motif for toe TTX/STX binding site of the ?*a+ 
channel— an andparaUel arrangement of ^-hairpins of repeats MV. Resi- 
dues that interact with toxins are marked by arrows. The numbered residues 
ore located in positions f + 1 of 0-tums. Relative orientation of 0-mms of 
the four hairpins is shown by flat arrows connecting the octal and bcta2 
strands, and it corresponds to an overhead view. 



binding cavity. In this model the residues at t 4- 6 can be 
ignored, because they are above the level of the toxin mol- 
ecule. This structural motif of the toxin binding site forms a 
funnel-tike aggregate barrel, consisting of four antiparallel 
P-hairpms wim the beta2 strands facing the vestibule and the 
betal strands on the outside toward the protein core (Fig. 4). 

Comparison of the structural motif to mutation 
data 

The arrangement of the four hairpins into an antiparallel 
structure results in a vestibule-like cavity appropriately sized 
for a toxin binding pocket. Although some of the mutation 



data were used to make our initial selection of interactive 
sites, model predicts correctly the results of mutation of all 
of the residues reported by Terlau et al. (1991). Neutraliza- 
tion of the aspartic and glutamic acid residues (but not swap- 
ping these residues) at positions i + 2 and i +> 5 of repeat 
I and i 4- 2 of repeat II resulted in loss of both TTX and STX 
block (Table 1). Neutralization of glutamic acid at position 
t + 5 of repeat II greatly reduces toxin block, similar to the 
effect of toxin analogs that affect the hydroxy! sites. The 
model also predicts possible participation of the residue at 
positions i 4- 3 in repeat I. Sarin et al. (1992) demonstrated 
that substitution of Phe 385 or Tyr in place of the naturally 
occurring Cys resulted in an increase in TTX and STX af- 
finity of about 1000 X. Backx er al. (1992) found that the 
reverse mutation of tyrosine to cysteine in the skeletal iso- 
fbnn reduced toxin affinity dramatically, and Hcinemann et 
aL (1992a) demonstrated the same effect for mutation of 
phenylalanine to cysteine in rBrIL Terlau et aL (1991) re- 
ported mutations in repeats I and II in positions i + 1 (Q383E, 
Q383K, G941E) and i + 6 (N388R) that had less than 10-fold 
effects on toxin block. The W386C mutation at £ + 4, which 
reduced TTX block 30-foid (Tomasclli et aL, 1993), is a very 
nonconsexvarXve change, and it undoubtedly caused some 
structural reorganization of the vestibule. Mutations beyond 
the lOamino acid 0-hairpins modelled here (R379Q, Q391K, 
> D949N, E952Q) had littfe effect (Terlau et al., 1991; Kontis 
and Goldin, 1993), suggesting that the dimensions of the 
proposed vestibule are appropriate. 

The mutants of residue Asp 1426 in repeat III and of residue 
Asp 1717 in repeat IV represent a good test of the antiparallel 
alignment of the 0-hairpins. This model predicts that side 
chains in position i + 5 of repeats HI and IV are directed 
outside the pocket. Indeed, mutants of this carboxyl of repeat 
ffl (D1426N, D1426Q, D1426K) do not affect toxin block 
significantly (Table 1). On the other hand, substitution by a 
positive charge at position i + 4 of repeat in did interfere 
with block of both toxins (M1425K). In repeat IV position 
/ + 4 is active, and all three mutants of Asp 1717 — D1717N, 
D1717Q, and D1717K do not bind STX. Besides, the residue 
at position i + 1 can influence toxin block (A1714E). In 
summary, the structural model we propose is consistent with 
the toxin sensitivity of the mutants of the SS1-SS2 peptides 
so fax reported. 

ENERGETIC IMPLICATIONS OF THE TOXIN 
BINDING POCKET 

Molecular modelling 

Modelling was accomplished in the Insight and Discover 
graphical environments (Biosym Technologies, Inc., San 
Diego). The molecular mechanics energetic calculations uti- 
lized the force field cvff (consistent valence force field), 
which was specifically derived to fit structures of peptides 
and proteins. In this approximation the formation of hydro- 
gen bonds is modelled as an electrostatic interaction. For 
minimi ration procedures the steepest descents and conjugate 
gradients have been used. 
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FIGWE 5 Spaxial mode] of the complex of TTX with the betal and beta2 strands of rcpeos I and n of the cardiac Na* channel fcofann. Four critical 
caiboxyl residues and the TTX arc ideaufied as space-filling imagea, while the remainder of the p-hairpins are shown with a backbone ribbon and stick 
sice chains- 



Modelling of the outer vestibule of the Na + channel was 
accomplished in three stages. First, the optimal conforma- 
tions of the ^-hairpins were calculated. Second* the inter- 
action of TTX with the repeat I and II peptides were obtained. 
The peptides of repeats HI and IV were assumed to form only 
the remaining walls of the pocket Finally, STX was sub- 
stituted into the repeat l-II structure derived from interaction 
with TTX and repeats III and IV were added. The optimal 
Spatial structures of TTX and STX were determined in the 
same program environment on the basis of their chemical 
structures and specific configurations (Kao, 1986). 

The reasons for modelling the decapepride of each repeat 
as a p-hairprn have been discussed, but each residue con- 
formation must be determined- In this structure the rwo 
/3-strands are united by a p-hira. As zero approximation for 
angles of rotation <KC" - N) and t^C" - C) in the backbone 
of the peptide chain of P-strands the angles of rotation in 
standard anriparallel ^structure, <f<C° - N) = -139* and 
,^C* - C) » +135° (Creighton, 1993) have been taken, 
while for the rwo central residues of the 0- turns with numbers 
i + 1 and i + 2 we used angles of rotation in reverse 0-rurns 
of type Hi: c^C* -N) = -60° and tfC - C) = -30* 
(Wilmot and Thornton, 1988). Angles <f> and if/ in 0- strands 
correspond to region B of conformational maps 4> - $ (<fr < 
0°, > 0*), while the angles (f> and *|r in 0- turns of type III 
correspond to region R (<f> < 0*. ift < 0°). Therefore we note 



the conesponding conformational states of the amino acid 
residues by the symbols & and R. 

As a preliminary step during the search for optimal con- 
formation of the 0-hairpins of repeats I and Q, minimization 
of potential energy of the pentapeptide fragments on the N 
and C ends in the conformations BBBBR and RBBBB was 
carried out. Since only the C ends of the rwo j3-hairpins (the 
beta2 strands) interact directly with TTX, we found the op- 
timal conformations of repeat I (Asp-Cys-Trp-Glu-Arg (for 
the cardiac channel), Asr>Tyr-Trp-Glu-Asn (for the muscle 
channel), and Asp-Fhe-Tip-Glu-Asn (for the rat brainU 
channel)) and of repeat II (GIu-Trp-De-Glu-Thr). These pen- 
tapeptides were allowed to interact with TTX and new ap- 
proximations of the angles of rotation x x (C* - C 0 ), ;&(C* - 
CO were obtained for the side chains of residues of Asp and 
Glu such that they were in maximal proximity to the gua- 
nidinium group of TTX. It was possible only when fragments 
Asp-X-X-Ghi and Gm-JT-Jf-Glu accepted the conformation 
RBBB; that is, the conformation of the beta2 strands. After 
minimization of the energies of these complexes the carboxyl 
groups of the side chains of Asp and Glu formed optimal 
hydrogen bonds with the guanidmium gronp of TTXL 

The optinial conformations of the BBBBR fragments of 
the N ends and the complex conformations derived above for 
RBBBB of the C ends of the 0-hairpins were used to form 
the two decapepdde 0-hairpins of repeats I and II in the 
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FIGURE 6 Spadal model of die complex of TTX with tie betal and beua strand* of repeals. I ud Q of the adult skeletal muscle isofonn of the Na*" 
channel. The residue Tyx 3 * 5 is also shown In a spacefilling image, Note die proximity of the aroma tie zing of Tyr to the hydrophobic surface of the TTX 



conformation BBBBRRBBBB. Fof angles of rotation $ and 
^ in residues i + 1 and t + 2 in the J^tums the initial con- 
dition that there was a hydrogen bond 0=0 — H — N be- 
tween residues i and i -h 3 in the chain was used. 

Finally, the three pair structure of 0-hairpins of repeats I 
and n and TTX was obtained by combining the TTX mol- 
ecule in its complexes with the ^-hairpins of both repeats 1 
and II for energy minimization of the ensemble. This re- 
sulting complex is shown for the cardiac isofonn in Fig. 5. 
The two /3-bairpins of repeats I and II form a cavity with the 
dimensions of TTX. The depth of the cavity is 12-13 A. The 
bottom is formed by two carboxyl groups of the side chains 
of Asp 3 * 4 and Glu* 42 in their i + 2 positions on the two 
/3-turas. The distance between the closest oxygens is 4.7 A, 
providing room for water molecules in between for com- 
pensation of electrostatic repulsion. The distance between 
the outer carboxyU of Glu 337 and. Glu* 45 are significantly 
larger at 8.1 A, so that they could not form direct hydrogen 
bonds with each other. The back wall of this cavity is formed 
by the side chain of Trp 943 of repeat n. As described before, 
three carboxyls simultaneously form salt bridges with 
the guanidinium group (Glu 387 with N(l>— H, Asp 384 with 
N(2) — H, and Glu* 12 with N(3>— H), while Glu 945 forms two 
optimal hydrogen bonds with the two hydroxyl groups 
C9-OH and C10-OH of TTX (Fig. 5). The distance between 
oxygens of the carboxyls and protons of the corresponding 
N — H groups of TTX in this complex are 2.4-2.6 A. These 
distances axe slightly greater than the corresponding dis- 



tances in complexes of TTX with individual chains of repeats 
I or II (2.1-2-2 A). In the three part complex the salt bridges 
have slightly less energy than with only one repeat and TTX. 
On the other hand, the distances O - • O in hydrogen bonds 
formed by side chains of Glu 945 with C9-0H and C10-OH 
groups are 3.0 A, corresponding to optimal hydrogen bonds 
OH — O ((O — O) = 2.8 ± 0.1 A (Pimentel and McCiellan, 
I960)). 

Calculation of the energy of interaction 

Nonbonded and electrostatic interactions were calculated. 
For the complex of repeats I and II and TTX for the cardiac 
isofonn, energies of nonbonded interactions with residues 
384, 385, 387, 942, and 945 (rBrll numbering) were -1.5 to 
-2 kcal/mol, although Trp 943 bad an energy of -4.8 kcal/mol. 
Electrostatic interactions were calculated using a dielectric 
constant e ~ 10, appropriate to calculations for peptide 
chains in polar (water) media (Krimm and Mark, 1968). En- 
ergies of electrostatic interactions for each of the glutamic 
and aspartic acids, including Glu 945 with two hydrogen 
bonds to C9-OH and C10-OH is about -7 kcal/mol. In ad- 
dition to these, there is electrostatic interaction with the posi- 
tive side chains of Arg (378, 388, and 937), although there 
are no nonbonded interactions with them.' Therefore, the total 
of nonbonded and elecrrostaric interactions of TTX with the 
0~hairpins of repeats I and II of the cardiac isofonn are cal- 
culated to be -15.6 and -18.2 kcal/mol, respectively. 



| 2003.11.13 09:10 ^0771+3833520 — — — — ZX* Q 



Biophysical Journal 



Volume 66 January 1994 



1 I . ... 



7 Spatial model of the complex of STX with the betal and beo*2 strand* of repeats MV of the cardiac Na" 
of the second guanidinium group and carbamyl of STX id the center of the model to Asp 1717 of repeat IV. 



."*" channel tsofomv Note the proximity 



The above calculations were made for the cardiac isofonn 
so that we could compare the model to the experimental 
results of Satin et al. (1992). They reported that substitution 
of a Tyr or Phe for Cys 385 increased the sensitivity of TTX 
block about 1000 times. This suggests that there is an im- 
portant nonpolar interaction with TTX. Indeed, TTX has 
asymmetrical distribution of -OH groups. The surface op- 
posite to tbe C9 and CIO hydroxyls contains mainly protons 
of C — H groups. Replacement of the Cys 385 with either Tyr 
or Phe causes the aromatic rings to be parallel to the nonpolar 
surface of TTX (Fig. 6). The nonbonded energy of interaction 
for tyrosine or phenylalanine is -7 kcal/mol, while that for 
cysteine is only -2.1 kcal/mol, and energetic gain of about 
5 kcal/mol for interaction with the repeat 1 0-hairpin, as sug- 
gested by Satin et al (1992), and generally supports the struc- 
tural motif. 

The 7,8^-guamduuum group of STX forms the same 
physical relationship and energy of interaction as that of tbe 
TTX guamdinium group. As before, this group takes part in 
the same salt bridges with Asp 384 , Glu 387 , and Glu* 42 - The 
interaction of Glu 945 with two hydroxyls is preserved by the 
two hydroxyls at C12 of STX. After energy minimization of 
the three element complex, optimal hydrogen bonds were 
formed with O ••• O distances of 3.0 A. It is apparent that 
STX shares the same sites as TTX on repeats I and II in 
accordance with the close agreement with their biological 
action (Kao, 1986). The nonbonded interaction energies are 
almost identical. However, the second 1,2^-guaiudinium 



group influences the electrostatic interactions in two ways. 
The energies of electrostatic interaction of STX with the side 
chains of Asp and Glu (384, 387, 942, and 945) are larger 
in absolute scale that those for TTX, but this is compensated 
by repulsion of the Arg interactions. Consequently, the en- 
ergies of binding of TTX and STX with the two ^-hairpins 
of repeats I and 12 are almost identical 

The other two 0-hairpins of repeats HI and IV were aligned 
by allowing a salt bridge between the 1,2,3-guamdimum 
group of STX and the side chain of Asp 1717 of repeat IV and 
the presence of van der Waals contacts between STX and the 
side chain of Met 1425 of repeat HI. Neutralization of Asp* 717 
(D1717N) reduced STX block. Substitution of Gin for 
Met 1425 had little effect, but its replacement by the positively 
charged Lys prevented block (Terlau et al, 1991). We in- 
terpret this to mean that the side chain of Met has a non- 
bonded interaction with STX Finally, the complex of STX 
and the four 0-hairpins of repeats 1-IV was evaluated, in- 
cluding dense packing of p-hairpins but omitting van der 
Waals repulsions of their neighboring side chains (Fig. 7). In 
this structure two protons at N(2) of the 1^-guanidinium 
group form hydrogen bonds with the carboxyl of Asp 1717 , 
and the carbamoyl group of STX also takes part in effective 
electrostatic interactions with this residue of repeat IV. Al- 
together, the energy of nonbonded and electrostatic interac- 
tions of STX with the fj-hairpin of repeat TV is -13 kcal/mol. 
The energy of interaction of STX with the ^-hairpin of repeat 
HI is significantly smaller -1.4 kcal/mol. The Cartesian co- 
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ordinates of some of the important groups of this complex are 
given in Table 3. 

The difference in energy of interaction between STX and 
TTX was estimated by substituting TTX for STX in the tour 
3-hairpia model. The nonbonded interaction energies with 
repeats III and IV are the same, but the loss of the additional 
electrostatic energy of the 1,2^-guanidinium group reduced 
the TTX interaction energy by 6 kcal/moU giving it a lower 
affinity than STX. We did not calculate free energy of these 
interactions* so these values cannot be directly converted into 
binding affinities for experimental comparison. However, 
the calculated energy differences between isoforms and be- 
tween toxins are experimentally consistent 

Although the proposed model of the Na"*" channel vesti- 
bule as the toxin binding pocket fits the mutation data pres- 



ently available, there are several uncertainties to be consid- 
ered. First, the outside of the pocket faces other parts of the 
protein, which could have significant effect on the shape and 
stability of the /3-hairpin region. Second, the activities of 
water and protons in the vestibule could be restricted, af- 
fecting the electrostatic interaction energies and the local pH. 
Finally, the effects of such mutations as K1422E and 
A1714E (Terlau et aL. 1991) do not mean that in native 
channels Lys and Ala interact with the toxins, because the 
addition of carboxyls at these sites may create new loci for 
toxin interaction. 



0-TURNS OF THE TOXIN BINDING POCKET . 
FORM THE ION SELECTIVITY FILTER 

The mutation studies of Terlau et al (X991) included esti- 
mates of the single channel conductances from noise analy- 
sis. Although many of the mutations had no effect on single 
channel conductance, neutralization of the glutamic acid or 
aspartic acid residues reduced iu Specifi c ally, neutralization 
of Asp 3 * 4 (D384N) and Glu* 42 (E942Q) greatly reduced 
ionic current (Table 1). These residues are located in posi- 
tions i + 2 of 3-tums of repeats t and II, and they axe ap- 
parently involved both in toxin binding and ion conduction 
of the Na* channel. Subsequent studies by Heinemann et al. 
(1992b) showed that simultaneous mutations of two residues 
in repeats m and IV (K1422E and A1714E) made the Na* 
channel permeable to Ca 2+ . They suggested that the resulting 
four negatively charged side chains at positions 384, 942, 
1422, and 1714 reproduce the structure of the selectivity 
filter of Ca 2 * channels. Such a change in selectivity is in 
accordance with our structural model, because the two mu- 
tations refer to residues at positions i + 1 of the p-turns of 
repeats HI and IV, which because of this antiparallel arrange- 
ment are directed to the inside of the binding cavity (Fig. 8). 
It is interesting to consider the hypothesis that two carboxyls 
at this site are responsible for Na* selectivity of the native 
Na* channel, but that three or four are required for Ca 2 * 
permeation. 

The ions Na* and Ca 2 * have the same ionic radii (1.02 and 
1.00 A) in crystallohydrates. They are typically octahedrally 
coordinated with six molecules of water, but because of the 
difference in charge, they have very different hydration en- 
ergies (-105 kcal/mole for Na* and -397 kcai/mole for 
Ca 2 *) (Hille, 1992). This leads to the idea that electrostatic 
interaction of Na* with the two carboxyls at positions 384 
and 942 are able to reduce the hydration sheU energy suf- 
ficient to decrease the barrier to Na* movement through the 
channel pore. However, three carboxyls may be required for 
weakening of the hydration shell of Ca 2 *. Three are sug- 
gested because the single mutant K1422E and the double 
mutant IC1422E and A1714E have the same Ca 2 * conduc- 
tance (Heinemann et aL, 1992), so that only the carboxyl at 
position 1422 may be involved. 

In the four /S-baiipin model of the Na* channel with the 
toxin removed, Na* coordinates with Asp 394 and Glu* 42 to 
force out two or more molecules of water. Coordination of 
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FIGURE 8 Double mutant of the toxin binding site of die Na + ctand structure with tnxm removed (Kl and AJ714E (Heinemenn et al., 1992b)). 



Na"*-4H 2 0 in this region, which weakens binding energy 
between Na * and water in the transition state is shown in Fig. 
9. For the double mutant K1422E and A1714E (Fig. S), the 
carboxyls of three side chains (Asp 38 *, GIu 9 ** and Glu 14ZZ ) 
can be coordinated with Ca 2 "*"-3H 2 0, substituting for three 
molecules of hydration water. In the double mutant model 
Glu 17 ** does not coordinate with Ca 2+ because of the asym- 
metrical structure of this region. This region of the model is 
similar in size to that suggested by Hille (1971) for the di- 
mensions of the Na* channel selectivity filter (3;1 X 5.1 A). 
Evidently, selectivity filters in Na* and Ca 2 * channels pro- 
mote the formation of ions with only one molecule of water 
(Na^-HsO and Ca 2+ H 2 0), which then pass through the 
channel pore. 

DISCUSSION 

The purpose of the model was twofold. Firstly, we wished 
to determine if the recently available mutation data from 
cloned Na + channels, in combination with the structures of 
the toxins and their analogs, would lead to a structurally 
satisfactory toxin binding site. Secondly, defining the site 
would characterize the structure of part of the permeation 
path, if the toxin binding site is in the path. Most investigators 
have assumed that the toxins block by physically occluding 
the channel, (Hille, 1975; Strichartz, 1984; Kao, 1986). How- 
ever, Green et aL (1987) have offered biophysical evidence 



chat the toxin binding site is at a distance from the permeation 
path. The mutation studies of Terlau et al. (1991), Satin et 
aL (1992), Backx et aL (1992), and Kontis and Goldin (1993) 
identified several amino acids in the SS1-SS2 regions of each 
of the four Na + channel repeats as critical for toxin binding. 
Further, Terlau et aL (1991) and Heinemann et al. (1992b) 
showed that mutations of several of these amino acids change 
single channel conductance and selectivity of the channeL 
This demonstration of structural overlap of the toxin binding 
site and the conductance/selectivity region is strong evidence 
chat the toxin site is indeed in the permeation path. 

The proposed model is supported by appropriate interac- 
tions with amino add residues identified by mutation studies 
and by consistency with the interaction of various toxin ana- 
logs. It provides a rationale for differences in binding be- 
tween TTX, which interacts directly with the SS1-SS2 seg- 
ments of repeats I and IL and that of STX, which additionally 
interacts with the segments of repeats tiff and IV. It also 
provides a basis for predicting the blocking affinity of new 
analogs and the effects of other mutations in the binding 
region. This is illustrated by the model prediction that the 
differences in toxin binding between the brain/skeletal and 
the cardiac isof orms is the result of a hydrophobic interaction 
with the aromatic rings of Fhe or Tyr, while Cys may be a 
spacer element that interacts only weakly. . 

An unexpected but welcomed outcome of the model was 
its prediction of the narrow region delineated by residues of 
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caxboxyto of Asp* 14 and Chi 5 ** of the cardiac Isoform. Na* is in violet The 



the 0-turns of the four 0-hairpins (384, 942, 1422* 1714). 
This region has been proposed as a candidate for the selec- 
tivity filter by Heinemann et al. (1992b) on tbe basis of mu- 
tation studies. When the toxin was removed from the model 
pocket formed by the four ^hairpins, the narrow region at 
the mm points had dimensions dose to those predicted by 
Hille (1971, 1992) on the basis of permeation of organic 
molecules of various sizes and shapes. Assuming that the 
molecule is stabilized in this conformation by the rest of the 
underlying protein, the carboxyl residues at positions 384 
and 942 on repeats I and II were located correctly for in- 
teraction with the hydrated Na + ion (Na-6H 2 0), displacing 
two water molecules and reducing the energy of interaction 
of the other waters of hydration. However, a third carboxyl 
on repeat IE, shown by Heinemann et al (1992b) by mutation 
sLudies to produce a channel with significant Ca 2+ perme- 
ation, was correctly Located in the putative selectivity region 
to provide an additional displacement of a third water mol- 
ecule in its hydration shell and decrease the energy of hy- 
dration of the other water molecules. We consequently sup- 
port the suggestion that this region may represent the 
channel'!* selectivity filter. Moreover, substitution of either 
Glu334 or Glul086 by Lys in the homologous selectivity 
positions in the human cardiac Ca 3 * channel a 1 -sub unit 
leads to a significant increase in Na* permeation through the 
Ca 2 + channel (Tang et al., 1993). Indeed, according to our 
model, the spatial arrangement of carboxyls of pairs of glu- 
tamic acids 334. and 677 or 677 and 1086, which are nec- 
essary for Na* permeation of the Ca 2 "* channel is approxi- 



mately the same as in the triangle of residues 334, 392, and 
1422 of the Na 4 " channel (Fig. 8). 

The appropriate location of the residues of repeat EH for 
simulation of a Ca 2 *** coordination site was the consequence 
of the antiparallel alignment of the ^-hairpins. This does not 
conflict with possible pseudo 4-fold symmetrical arrange- 
ment of the transmembrane segments S1-S6 of repeats I-IV, 
because the lengths of the segments connec ting beta2 to S6 
are different. In the case of repeats I and n, they are 9 residues 
each, whereas for repeats III and IV, they axe 19 and 30, 
respectively, and consequently could accommodate addi- 
tional loops in repeats m and IV. 

Another outcome of the model is the prediction of the 
structure of the outer vestibule of the channel, d efin ed as the 
region large enough to accommodate the bulky toxin mol- 
ecules. The size of this putative vestibule and the location of 
charged residues are consistent with the modelling studies of 
Dani (1986) and Cai and Jordan (1990). They provide spe- 
cific dimensions for future modelling and predictions of the 
outcomes of mutation studies on surface charge phenomena. 
This model also allows consideration of interactions between 
the channel and peptide toxins such as ^-conotoxin, which 
competes with TTX and STX in blockade of the adult skeletal 
muscle and eel electroplax channel isoforms (Becker et al„ 
1992). 

In this model the peptide fragments betal and beta2 of the 
P region create only the outer vestibule of the Na* channel, 
with the selectivity filler at its bottom. These peptides do not 
form the remainder of the transmembrane pore, in contrast 
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10 some previous models of a + and channels. A simple 
examination of the P regions shows that the amino acid com- 
position is very different, with many charged residues in the 
Na channel mouth. Consequently, it is likely that there will 
be significant differences in the structures of the pore regions 
of the two channel types. 

An important but unresolved question is the possible pres- 
ence of an electric Geld in this vestibule and its effect on toxin 
binding. Schild and Moczydlowski . (1991), Satin et aL 
(1992), Back* et aj. (1992), and Doyle et aL (1993) identify 
a site for divalent ion block within the electric field. Because 
divalent ion binding is competitive with toxin binding, this 
implies that there is a field drop within this structural element 
and predicts that toxin binding should show some voltage 
dependence. 

Although the proposed model is effective in its prediction 
of tsoform/mutation and toxin analog studies, the interface 
of these four 0-hairpins with the rest of the protein could not 
be considered It should be clear that such interactions 
would be important for the shape and stability of the ves- 
tibule. The model does provide a basis for selecting other 
mutations of these peptide segments, which could result in 
structural refinement, particularly with regard to the 
mechanism of selectivity. 
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Actions of Chiriquitoxin on Frog Skeletal 
Muscle Fibers and Implications for the 
Tetrodotoxin/Saxitoxin Receptor 

L. Yanc and C. Y. Kao 

From the Department of Pharmacology, State University of New York Downstate Medical 
Center, Brooklyn, New York 11203 

ABSTRACT Chiriquitoxin (CqTX) from the Costa Rican frog Atelopus chiriqutnsis 
differs from tetrodotoxin (TTX) only in that a glycine residue replaces a methylene 
hydrogen of the C-l 1 hydroxymethyl function. On the voltage-clamped frog skeletal 
muscle fiber, in addition to blocking the sodium channel and unrelated to such an 
action, CqTX also slows the activation of the fast potassium current in - 40% of the 
muscle fiber population. At pH 7.25, CqTX is as potent as TTX in blocking the 
sodium channel with an ED 50 of 3.8 nM. Its ED 50 's at pH 6.50 and 8.25 are 6 8 and 
2.3 nM, contrasted with 3.8 and 4.3 nM for TTX. These differences are attributable 
to changes in the chemical states in the glycine residue. The equipotency of CqTX 
with TTX at pH 7.25 is explainable by an intramolecular salt bridge between the 
amino and carboxyl groups of the glycine function, all other surface groups in TTX 
and CqTX being the same. From available information on these groups and those 
in saxitoxin (STX), the TTX/STX binding site is deduced to be in a pocket 9.5 A 
wide, 6 A high, and 5 A deep. The glycine residue of CqTX probably projects out of 
the entrance to this pocket. Such a view of the binding site could also account for 
the actions of STX analogues, including the C-ll sulfated gonyautoxins and the 
21-sulfocarbamoyl analogues. In the gonyautoxins the sulfate groups are equiva- 
lently placed as the glycine in CqTX, whereas in the sulfocarbamoyl toxins the 
sulfate groups extend the carbamoyl side-chain, leading to steric hinderance to 
productive binding. 

I NTRODUCTION 

Chiriquitoxin (CqTX) is an analogue of tetrodotoxin (TTX) found in the skin and 
eggs of Atelopus chiriquensis, a harlequin frog of the central highlands of Costa Rica 
(Kim, Brown, Mosher, and Fuhrman. 1975; Pavelka, Kim, and Mosher, 1977). In 
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those earlier studies, CqTX was known to retain the basic structure of TTX but to 
differ in the C- 1 1 position by having a large substituent. In mouse lethality assays 
(Kim et al., 1975) and on the isolated frog skeletal muscle fiber it had the same 
specificity of action and potency as TTX (Kao, Yeoh, Goldfinger, Fuhrman, and 
Mosher, 1981). However, insufficient material was then available for a full clarifica- 
tion of the chemical structure or a definitive study of its actions on specific ionic 
channels. Repeated attempts to collect more material were unsuccessful until June 
1988. Using current separation methods and a high pressure liquid chromatography 
TTX analyzer, CqTX was isolated and purified. Its structure has also been deter- 
mined (Yotsu, Yasumoto, Kim, Naoki, and Kao, 1990). 

We report here our studies of the actions of this new batch of CqTX on the 
functions of the sodium and potassium channels of the voltage-clamped, isolated frog 
skeletal muscle fibers. The significance of these studies extends beyond a mere 
comparison of CqTX with TTX. In recent years the structure-activity relations of 
some analogues of both TTX and saxitoxin (STX) have been clarified to a point 
where most active surface groups of these toxin molecules have been identified (for 
summary, see Kao, 1986; also Hu and Kao, 1991; Yang, Kao, and Yasumoto, 1992a; 
Yang, Kao, and Oshima, 1992&). Because the actions of CqTX are unique among 
TTX analogues, knowledge of its structure and actions permits us to draw some 
inferences on the probable physical dimensions of the TTX/STX binding site, and 
also to speculate on its possible location. 

MATERIALS AND METHODS 

CqTX was isolated from the skin of A. chiriquensis and its structure was determined by various 
spectroscopic methods, including mass spectrometry, 'H- and 1S C-NMR spectrometry, and 
infrared spectrometry. Some details in stereochemical configurations were confirmed by 
derivatization experiments. The chemical studies can be found in Yotsu et al. (1990). The 
structure of CqTX is shown in Fig. 1. In comparison with that of TOC, the only difference is 
that a methylene hydrogen of ihe C-1 1 hydroxymethyl group of TTX is replaced by a glycine 
function in CqTX, Stated differently, CqTX has all the structural features of TTX, except that it 
has an additional -CH(NH 2 )-COOH attached to the C-ll position. 

The biological studies were performed on short segments of single isolated frog skeletal 
muscle fibers in a Vaseline gap voltage-clamp method (Hille and Campbell, 1976). Details of 
the procedures and methods of data acquisition and analysis can be found in other papers (Hu 
and Kao, 1991; Yang et al., 1992a). The actions of CqTX on both / Na and / K were studied. 

Muscle fibers were held at -90 mV and depolarized by step increments until -20 mV 
beyond the £ Na . To reduce errors attributable to series resistance, experiments were conducted 
in a bath solution containing only 44 mM Na + . As a result, the £ Na was at +20 to +30 mV. The 
lower external Na* concentration does not alter the ED 30 of the toxin (Hu and Kao, 1991). 
Current traces showing maximum / Na and at £ Na were selected for detailed comparisons. The 
current at £ Na was I K . and was used in previous studies of other analogues of TTX and STX 
primarily to assess the health of the fiber. However, CqTX also slows the activation of the fast 
/ K . Thus, the current records taken at £ Na were also useful for determining the rate of activation 
of / K under various conditions. Time constants of the activation of the fast / K were determined 
on the same fiber in the control and toxin-affected states (see Results). / tests for significance of 
differences of these data were based on pair comparisons. 

The maximum / Na 's in toxin-affected conditions were compared with those in control (no 
toxin) states, and the normalized relation (7 Wn») was used for dose-response curves. The 
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small amount of 7 K present at this voltage docs not affect the ED 50 determinations (Yang et at., 
1992a). Least-squares linear regression lines were fitted to data in Hill plots of log (l-P)/P 
C - / WAnJ vs. log toxin concentration. Standard errors of estimate were also obtained. 
ED 50 's and their standard errors were taken at (l-P)f P = I (Hu and Kao, 1991). / tests for 
significance of differences in the sodium channel blocking effect were based on unpaired group 
comparisons. 

RESULTS 
General Description 

Although it is generally held thatTTX affects only the sodium channel, recent studies 
on two analogues of TTX modified in the C-6 position indicate that there may be 



OH 




TTX CH 2 OH 

CqTX CH(0H)CH(NH 2 )C00H 

FIGURE 1 Chemical structure of CqTX compared with that of TTX. Note that CqTX contains 
all the features of TTX except that it has an additional glycine moiety on C-l 1. In ITX, and 
probably also in CqTX, the C-10 -OH group has a pKa of 8.8. In CqTX, the C-12 amino 
group has a pKa of -9 and the C-l 3 carboxyl group has a pKa of ~2. These properties 
probably influence the sodium channel blocking action of CqTX at various pH's. 



some exceptions (Yang et al., 1992a). In about half of the muscle fibers from Rana 
temporaria, TTX and 11-deoxyTTX, which retains the general configuration of the 
groups around C-6, significantly slow the activation of the fast / K . 6-epiTTX, in which 
the -OH and the hydroxymethyl group of C-l 1 are in an epimeric configuration from 
that in TTX, does not cause significant slowing of the / K . Fig. 1 shows that CqTX 
shares the general structural features of TTX at the C-6 position; Fig. 2 shows that in 
some fibers it also shares the effect of slowing the fast / K . In such fibers, while / Na is 
markedly reduced and eventually blocked, the fast 7 K is also appreciably slowed. At 10 
ms of a depolarizing step the / K can be reduced by variable amounts, but at 20 ms it 
has usually reached the same amplitude as in the control state. 
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Effect of CqTX on Potassium Current 

Because 1 K is rather variable from fiber to fiber and even in the same fiber at different 
times, we used an arbitrary criterion to determine whether a toxin slowed the fast 7 K . 
The time constants (t) of the first 2 ms of the / K of the pretoxin and recovered states 
were separately determined and then averaged. The averaged control t was com- 
pared with that of the toxin-affected state. If the latter exceeded the former by 10%, 
we selected that fiber as having had its 7 K slowed. Table I gives a summary of these 
comparisons. Of 36 fibers used for this study, 14 showed some slowing of the 7 K . At 1 




Figure 2. Effect of CqTX on # Ka and 7 K of frog skeletal muscle fiber. Two different fibers (A 
and B from one; C and D from the other) in bath medium containing 44 mM Na + . Holding 
potential -90 mV. Currents produced by two voltage steps are shown, one showing maximum 
inward /, Va and the other, at £ Na , showing / K . In all panels, traces marked / denote the control 
state. In A and B, traces marked 2 denote responses to 4 nM CqTX, and traces marked 5 
denote responses to 8 nM CqTX. In C and D, traces marked 2 denote responses to 2 nM 
CqTX. and traces marked 3 denote the recovered state. A and C are at fast sweep; total 
duration is 10 ms. B and D are at slower sweep; total duration is 20 ms. In addition to reducing 
/.*., CqTX can also slow fast 7 K without appreciably reducing total / K . This feature is best shown 
in C and D, where the seemingly reduced / K at 10 ms (C) is no different from the control state 
at 20 ms (Z>). The effect of CqTX is nearly hilly reversible (C and D). In A and B, increasing 
concentrauons of CqTX increased reduction of / Na , but had no similar effect on the slowing of 
/ K . See text for additional reasons why the slowing of / K is independent of block of 7 Nj . 



nM CqTX generally did not cause any significant slowing of the 7 K , even though it 
could reduce / Na by -20% (see Fig. 5). In 2-40-nM concentrations CqTX often 
caused some slowing, but the degree of slowing was not dependent on the concentra- 
tion (Fig. 2, A and 5). In those fibers in which slowing occurred, 7 K usually recovered 
when CqTX was washed out. On subsequent application of another dose of CqTX, 
usually of a higher concentration, a lesser degree of slowing of 7 K often occurred 

No influence ofTTX. At pH 7.25, the ED 3( /s of CqTX and TTX are the same, 
and like TTX, CqTX does not affect the onset or decay of the 7 Na (see below). These 
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TABLE 1 

Effect of CqTX on the Rate of Activation of the Fast Potassium Current 



Time constant (t) of fast I K 

Concentration Fibers p 

Control Toxin 



nM 






ms 




\ 


12 


0.61 ± 0.02 


0.62 ± 0.02 


>0.05 


2 


10 


0.61 ± 0.02 


0.69 ± 0.04 


<0.05 


4 


22 


0.60 ±0.01 


0.67 ± 0.02 


<0.05 


8 


B 


0.63 ± 0.02 


0.73 ± 0.02 


<0.05 


16 


8 


0.45 £ 0.02 


0.56 ± 0.04 


<0.05 



All values are means ± SEM, with Fibers referring to the number of fibers used. 
Altogether 36 fibers were used, some of them for multiple doses of CqTX. Most fibers 
used for 16 nM came from frogs in May; fibers for all others came from frogs in 
July-September. P denotes P value in paired / test. 



observations suggest that the affinity of the two toxins for the binding site are similar. 
The following experiments on the interaction between CqTX and TTX were 
conducted with concentrations of TTX 1 0-100-fold higher than that of CqTX to skew 
the mass-action equilibrium in favor of TTX, Even under these conditions, the effect 
of CqTX on 7 K was not influenced by the presence of TTX, whether TTX was applied 
before or after CqTX. 

Fig. 3 is an example of a high concentration of TTX that did not prevent CqTX 
from slowing the fast 7 K . The 7 Na was fully blocked by 40 nM of TTX, and the fast 7 K 



1 2 3 4 5 

40% Na Ringer TTX 40 nM ™4J «|J TTX 40 nM 40% Na Ringer 




0.75 I 2 ms 
M 1 




Figure 3. Prior action of TTX has no effect on the slowing of 7 K by CqTX. (Top row) Current 
traces in various conditions sequentially as marked on top. Holding potential -90 mV; currents 
result from step voltages of -65, -40, -15, -10, and +35 mV. (Bottom row) Superimposed 
current tracings show different effects, with the number of the trace referring to the conditions 
marked above the top row. 40 nM TTX (panel 2) fully blocked / Na hut only slightly and 
insignificantly affected fast 7 K . An additional CqTX (panel 3) markedly slowed fast / K . This 
effect was reversible (panel 4). When TTX was removed (panel 5\ both 7 Na and / K reversed 
fully. 
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was slightly slowed (panels 7 and 2). When 4 nM of CqTX was applied, no additional 
eflects on 7 Na could be detected. The fast 7 K was appreciably slowed, and the 7 K at 10 
ms was substantially lower than before (panels 2 and J). When CqTX was removed 
the rate of activation of / K recovered to that in TTX (panel 4\ compare with panel 2) 
When TDC was also removed, both 7 Na and 7 K fully recovered (panel 5; compare with 
panel 7). This absence of influence of a prior application of TTX was demonstrated 
just as clearly in another fiber. Because the fibers were in very good condition, and 
the effects were clear, we did not repeat this sequence on more fibers. 

High concentrations of TTX did not reverse the effect of CqTX on the 7 K . In 10 
fibers, 7 Na was reduced by an ED 50 dose (4 nM) of CqTX. In every fiber the fast 7 K was 
appreciably slowed and the 7 K at 10 ms was reduced. Then a high concentration of 
TTX (100 nM in four fibers and 400 nM in six fibers) was applied. While 7 Na now 
became fully blocked, no further change in the 7 K was seen. In every case the toxins 
were allowed to act for at least 8 min to reach steady state. Since the effect of TTX on 
7 K was usually less than that of CqTX, if the high concentration of TTX had reversed 
the effects of CqTX on 7 K . one might have expected some recovery. Instead, the 7 K at 
10 ms in CqTX was 1.07 ± 0.11 (mean ± SEM of 10 fibers), compared with 
0.99 ± 0.1 1 uA in a combination of CqTX and TTX. 

No influence of [Na+] 0 or of direction of I Na . In addition to the apparently 
independent actions of TTX and CqTX on the 7 K , other evidence shows that the 
effect of CqTX on the 7 K is unrelated to functions of the sodium channel. In a 
Na*-free bathing solution, CqTX caused as much slowing of the 7 K as it did in 44 mM 
Na*. Moreover, slowing of 7 K was seen whether the 7 Na was inward or outward. 

The slowing of the fast 7 K by CqTX is qualitatively similar to that seen with TTX 
and 11-deoxyTTX in about half of the muscle fiber population (Yang et al., 1992a). 
From all the above observations, we conclude that the slowing of the fast 7 K is not 
related to the sodium channel blocking action of these toxin molecules, but that it is 
somehow related to the configurations of the C-6 groups common to all three. See 
Discussion for further comments. 

Effect of CqTX on Sodium Current 

Similar mechanism to that of TTX. Because CqTX slows the 7 K in an appreciable fraction 
of the fiber population, a question arises as to whether it might have slowed the 
mactivation of 7 Na such that some overlap of currents might be responsible for 
slowing the 7 K . To examine this problem, experiments were performed on fibers in 
which the 7 K had been blocked by applying Cs+ to the end pools and allowing it to 
diffuse into the nodal area. Fig. 4 shows a typical example of such an experiment. In 
the control state, only the inward 7 Na was seen. Upon action of 4 nM of CqTX the 7 Na 
was reduced but there was no prolongation of the decay phase. When the toxin- 
affected 7 Na is scaled to match the peak of the control 7 Na , the currents are 
superimposed, with no detectable differences in either the onset or decay kinetics 
I he same results were seen in nine other fibers. Thus, the mechanism of action of 
CqTX on 7 Na is apparently similar to that of TTX. 

Dose-response relation. Fig. 5 shows the dose-response relation of CqTX on the 
/ Na at pH 7.25. As in the case of TTX, the data fit well a bimolecular reaction scheme 
where one toxin molecule interacted with one receptor site. This point is reinforced 
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Figure 4. CqTX docs not affect the rate of inaciivation of / Na . Two different fibers, both with 
/ K blocked by Cs\ Holding potential -90 mV; step voltages to +10 mV. Traces marked / are 
maximum / Na in control state; traces marked 2 are maximum / Na in 4 nM CqTX (EO 50 ). (fnsels) 
Traces 2 up-scaled to have peaks coincide with those of traces /, showing the absence of effect 
of CqTX on the kinetics of onset or decay of 7 Nil . 



by the slope of 1.1 in the Hill plot where log (\-P)fP (P = 7' Na // Na ) is plotted against 
log toxin concentration. The ED 50 of CqTX of 3.8 nM is not statistically significantly 
different from that of TTX, which is 4.1 nM (see Yang et a!., 1992a). 

Effect ofpH on the dose-response relation. Because various toxin molecules contain 
ionizable groups, alterations in pH have been useful in correlating chemical changes 
with potency differences. Also, because the actions of CqTX are rather similar to 
those of TTX, it is instructive to examine the influence of pH on the actions of CqTX 
in comparison with those of TTX. In TTX, the C-10 -OH has a pKa of 8.8. It is 
more potent at pH 7.25 than at pH 8.25, with a relative potency that agrees well with 
the relative abundance of the protonated species of the C-10 hydroxyl group (Hu and 
Kao, 1991). Although there is no direct information on the property of the C-10 




Figure 5. Dose-response relation of CqTX in blocking / Na at pH 7.25. {A) Relative residual 
' Na vs. log concentration of CqTX. The solid curve represents a bimolecular reaction scheme. 
The number attached to each point denotes the number of individual fibers from which data 
are obtained. Symbols denote mean values and vertical bars denote I SEM if larger than the 
symbol. The ED 50 of 3.8 nM is not statistically different from that of 4. 1 nM TTX (see text). (B) 
Hill plot, where P = /' Na // N v EO 50 is determined when (1-P)/P = 1 . A slope of 1 . 1 is consistent 
with a bimolecular reaction. 
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group in CqTX. because of similarities in chemical structures to TTX one might 
assume the pKa to be about the same as that of TTX. On this basis, one might expect 
CqTX to be similarly more potent at pH 7.25 than at 8.25. 

Fig. 6 shows that, in actuality, CqTX is appreciably more potent at pH 8.25 (ED^, 
2.3 nM) than it is at pH 7.25 (ED 30 , 3.8 nM; P = 0.01). Also, at pH 8.25 CqTX is 
significantly more potent than TTX (ED 50 , 4.3 nM; P < 0.001, group comparison 
with data in Hu and Kao, 1991). These differences are probably not due to 
differences in protonation of the C-10 group, but are attributable to dissociable 
groups in the glycine function of CqTX. 

At pH 6.50, increments in the protonated species of the C-10 group over that 
present at pH 7.25 are small, and the potency of TTX is slighdy lower (ED 50 , 3.8 
nM), possibly because of a decrease in the density of surface negative charges (see 
discussion in Hu and Kao, 1991). For CqTX, however, the potency at pH 6.50 (ED 50 , 




Figure 6. Influence of pH on the dose-response relation of CqTX in blocking / Na . Conven- 
tions are similar to those in Fig. 5, from which data at pH 7.25 are taken. Slopes in the Hill plot 
are 0.9, 1.1, and 1.1, and ED 50 's are 6.8, 3.8, and 2.3 nM for pH 6.50, 7.25, and 8.25, 
respectively. See text for explanations of different potencies. 



6.8 nM) is significantly less than that at pH 7.25 (ED 50 , 3.8 nM; P = 0.05). Again, the 
reason for the difference must be sought in changes in the states of the glycine 
function. 

DISCUSSION 

Until the discovery of CqTX in 1975, no analogue or derivative of TTX was known to 
have anything but a small fraction of the potency of TTX. Thus, the chemical 
mechanism by which TTX blocked the sodium channel remained largely in the realm 
of speculation. CqTX was the first analogue of TTX to have substantial potency, and 
it rekindled studies of the structure-activity relations of the TTX family of molecules. 
Coincidental^, natural analogues of STX were also being discovered. The availability 
of analogues and derivatives of both TTX and STX fostered a reexamination of the 
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chemical mechanism of the sodium channel blockade. As shall become clear pres- 
ently, because of its structure and its actions CqTX is of special importance in helping 
to consolidate a conception of the probable physical dimensions of the TTX/STX 
binding site. 

Relation to Previous Work and Effect on I K 

The first study of the cellular actions of CqTX was made at a time when no 
voltage-clamp method capable of controlling the rapid sodium current was available 
for work on the frog skeletal muscle fiber (Kao and Yeoh, 1977). The maximum 
dV/dt of the action potential was used for an indirect assessment of sodium channel 
function, and constant current methods were used to study potassium channel 
function. In those experiments on muscle fibers of Rana pipie-ns, slowed repolarization 
of the action potential and a reduced level of delayed rectification were seen in some 
fibers at the same time that spike generation was impaired (Kao and Yeoh, 1977; Kao 
et al., 1981). Based on changes in the maximum dV/dt of the action potential, CqTX 
was judged to be as potent as I I X in blocking the sodium channel, but CqTX was 
also thought to have some effects on the potassium channel. The present results 
confirm that CqTX does have some effect on the potassium channel, but show that 
the effect is not a simple blockage as first believed. It is a more complex interference 
with the activation kinetics of the fast / K and occurs in only about half of the 
population of muscle fibers. Because of variability, we do not have enough informa- 
tion at this time to fully explain that effect. All evidence in hand indicates that the 
effect on is unrelated to the sodium channel blocking action. The effect on /k is 
also seen with TTX and 1 1-deoxyTTX, but not with 6-epiTTX or neoSTX (Yang et 
al., 1992a). The common structural feature shared by CqTX, TTX, and 11- 
deoxyTTX is the configuration of the groups on C-6. Perhaps some potassium 
channel molecules, possessing some homology with the sodium channel (eg., Tem- 
pcl, Papazian, Schwarz, Jan, and Jan, 1987) have a binding site for those groups. 

Effect on l^a and Influence of pH 

In the sodium channel blocking effect, CqTX is unusual among TTX analogues in 
possessing a high degree of potency. With the exception of 1 1-oxoTTX (Wu, Yang, 
Kao, Yotsu, and Yasumoto, 1991), all other TTX analogues studied have varying 
degrees of reduced potency as compared with that of TTX, but CqTX is equally as 
potent as TTX. However, the pH influence on the potency of CqTX is rather 
different from that of TTX. Based on known chemical properties of CqTX (Fig. 1) we 
offer the following explanations for these differences: 

In the glycine moiety of CqTX, the C-12 amino group has a pKa of ~9, and the 
terminal (C-13) carboxyl group has a pKa of -2 (Yotsu et al., 1990). At pH 7.25 the 
C-12 amino group is mostly protonated, whereas the C-13 carboxyl group is 
deprotonated. An intramolecular salt bridge between those vicinal groups could 
largely remove them from any charge interactions with receptor groups in the 
sodium channel. As all the other groups in CqTX are the same as those in TTX, 
including the C-l 1 -OH group that forms a hydrogen bond with a receptor site (Yang 
et al., 1992a), the equipotency of CqTX and TTX is not difficult to understand, 
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provided a satisfactory explanation can be found for the absence of steric influences 
of the glycine moiety (see below). influences 

Tl^and ^i"""^ T P ° tenCy S ° different between THC and CoTX> For 

X^T? H g yC J nC fUn ?° n h3S "° addkional ^on-able groupTsen^ ' to 
pH changes. Under acidic conditions, the C-13 carboxvl sroun in r„TY 7 T 
form an intramolecular lactone with the C-6 -OH « ^990 22 r n 

^though pH 6.50 is less drastic than the co^2ZtZ^Z l c 
may be a change in the equilibrium between CaTX and PoTY u r 

b" h C 6 and c 1 1 OH P ^ " * dilfcra " P"* 0 " *»" » TTX. As 

Implications for the TrX/STX Receptor 
In addition to electronic factors, steric factors also influence the docking orientation 

~o analogue* nxh of TTX (6*,,TTX J niS^TTO Z fjr^u' "T h 

^tfSTr Carbam0yl fUnCti0 " ° f STX ^ponded ,0 th™! 1 

crevice or a fold of the uJL^SJ^^ ^ "* ™» be ^ d * - 
Combmrng the present results with the earlier deduction, the TTX/STX binding 
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site can be viewed as being a pocket ~ 9.5 A wide, 6 A high, and 5 A deep (Fig. 7; also 
Kao and Yang, 1992). In it are seven site points (designated sites d-g), which interact 
with various surface groups of either the TTX or the STX (and neoSTX) molecule. 
The CqTX molecule could have an additional reactive site point complementary to 
the C-12 -NH 2 . However, the most impressive feature of CqTX is that in spite of the 
additional steric volume provided by the glycine function, it is very similar to TTX in 
potency. One possible explanation is that when bound productively to the receptor, 
the CqTX molecule is so oriented that the glycine function faces, and possibly 
protrudes out of, the entrance to the binding site, hence the minor contributions of 
the additional reactive groups on CqTX. 

An attractive aspect of such a view of the binding site is that it can also explain two 
features regarding STX analogues that have not been satisfactorily accounted for 
before. Since TTX and STX bind to the same receptor, this ability to accommodate a 
group of chemically different compounds adds to the general plausibility of the 
proposed view. Among STX analogues there arc two sulfur-containing subfamilies, 
the C-ll (numbering system for STX) gonyautoxins, which contain large sulfate 
substituents on the C-l 1 position, and the 2 1 -sulfocarbomoyl toxins, which contain a 
sulfate end on the carbomoyl side-chain (for summary, see Oshima, Sugino, and 
Yasumoto, 1989). For gonyautoxins, the -OSO* substituents on C-l 1 add both large 
steric volumes and strongly negative charges rather close to the critical C-12 
gem-diols, yet their potencies are only minimally reduced (Kao, Kao, James-Kracke, 
Koehn, Wichmann, and Schnoes, 1985). On the other hand, the bulk and charges of 
the -OS0 3 groups in the sulfocarbomyl toxins at the end of a side-chain profoundly 
lower the potency. Both of these features can be explained by the proposed view of 
the TTX /STX receptor as situated in a pocket. 

Molecular modeling shows that when the active guanidinium groups of TfX and 
STX are aligned, as are the C-9 and C-10 -OHs of TTX with the C-12 -OHs of 
STX, the C-ll sulfate groups of the gonyautoxins are stereochemical^ placed much 
like the glycine moiety of CqTX (Fig. 7). Perhaps their minimal influence on potency 
is also to be attributed to their projecting out of the entrance to the binding site 
pocket and being removed from nearby reactive site points. On the other hand, the 
sulfate groups on the 21-sulfocarbamoyl toxins extend the carbamoyl side-chain so 
much that steric hinderance interferes with the entry of the toxin into the binding site 
pocket. Understandably, when those sulfate groups are cleaved through acid hydrol- 
ysis (Shimizu, Kobayashi, Genenah, and Oshima, 1984), the product toxins become 
considerably more active because they are now freed of the steric interference and 
can productively occupy the binding site to block the sodium channel. 

Observations on mutated sodium channels of rat brain suggest that negatively 
charged amino acid residues in short segments 2 between S5 and S6, at positions 384 
and 387 of repeat I, and equivalently placed residues in repeats II, III, and IV 
strongly influence the actions of TTX and STX (Noda, Suzuki, Numa, and Stuhmer, 
1989; Kontis and Goldin, 1991; Terlau, Heinemann, Stuhmer, Pusch, Conti, Imoto, 
and Numa, 1991). Mutating the anionic glutamate or aspartate to the uncharged 
glutamine or asparagine raises the ED50 for blocking / Na by more than three orders of 
magnitude. We speculated that if glutamate 387 corresponded to our binding site a, 
then the carbonyl oxygen of asparagine 388 could represent our binding sites b and c 
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(Yang et al., 1992a, b). However, the situation is complicated by the contradictory 
observations that mutants containing arginine 388 showed no significant difference 
in sensitivity to TTX and STX (Terlau et al., 1991), whereas in the sodium channel of 
rat cardiac myocyte TTX/STX insensitivity is associated with the presence of arginine 
at 388 (Rogart, Cribbs, Muglia, Kephart, and Kaiser, 1989). Not unexpectedly, 
molecular models of short segments of amino acid residues from 383 to 390 do not 
show any folds that could envelop the TTX or the STX molecule. However, in the 
concentric model proposed by Guy and Conti (1991), the TTX/STX binding site 
could straddle two repeats, which could contain the appropriate residues to provide 
the crevice that serves as the toxin binding site. Experimental identification of the 
binding site would require specific marker compounds that need to be made. 
However, a successful identification of the binding site could add much to a better 
understanding of the structure and function of the sodium channel. 
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The Guard dinivtm Toxin Binding Site 
on the Sodium Channel 
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Summary 

Study of TTX and STX toxins and their interacaon with cloned Na 
channels has led to a molecular model of the toxins' binding site. This model is 
able to explain isoform differences in binding affinity, allowing prediction ^f the 
structures or toxin-resistant channels. It provides an. example of a detailed drug 
binding site that could serve as an example for drug engineering to improve 
affinity and specificity. The model has also suggested important molecular 
characteristics of the channel's permeauon path, providing a wealth of ideas lor 
study of this fundamental biophysical process at the molecular level. The model 
of the TTX and STX binding site outlines an approach to resolution of mo- 
lecular structure :.nd of certain structure-function reladonships, and "V™"}" 
tation provoked by the model will assist in improving it. (Jpn Heart J 1996; 37: 
(i»3-(>92) 

Keywords: Guanidinium toxin Tetrodotoxin Saxitoxin Sodium 
channel Molecular modelling Vestibule Structure-function 

THE Ruunidinium toxins, tetrodotoxin (TTX) and saxitoxin (STX), are 
highly poisonous marine toxins that have become extraordinarily valuable 
tools in the study of excitable tissues. They bind uniquely and with high affinity 
to voltage-gated sodium channels and block them completely. The usual vector 
for poisoning for TTX \. ^Bj^SkLMh and brJ&L^s™^ -th 
Ingestion of cith^xin-p^ HowcVcr ' ?* f^. 

does not s^hc^XDi^ instcad_CQncentratcs..it via l*acte.ria in its food, 
cESnTStarijSS^^ are mgestcd by the 

s^lfisTO^hcr animals can also concentrate the toxins, presumably derived from 

their food. m m , n . . m . _ 

Biophysical evidence points to a stoichiometric binding of the toxins to a 
shared site in the outer mouth or vestibule of the sodium channel pore and party 
within the membrane field. Mutational studies with cloned mammalian sodium 
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channels has identified an involved str^K „r ... 

protein domains of the a-subunh It h" h ° f amm ? I aads ,n each * ** ^ur 

model of this bindin.site 1 P ° SS,bIC t0 a m ° IeCular 

rerion resu tin* t I T Wh,Ch appcars to ^iude the channel's selectivity filter 

Xelec ^7 " PrCdlCtl ° nS ab ° Ut ^ m ° lCtU,ar of ion pe Jeaticl 

The GuANmiNruM Toxins 

299 ^ S Y 3 u Sma " rfgid ^ m ° ,CCU,es CTOC = 319 MW; STX = 

299 MW). They have been synthesized and structurally determined by X rav 

ToZ in", h" C T mhCdCal,y ' fadHtat ^ id -«Sadon of cnucal acdve 
groups mvolved m toxm binding.- For TTX the most critical groups are Ae 
,2 3 guan.duuum and the C9 and CIO hydrops, and for STX they a c die 
7,8,9 guamdmxum and the two G12 hydroxyls. Presumably the 7 2 3 
guamd.mum .s also important, but the analog evidence is less clear" W 
Narahash, and colleagues" demonstrated that TTX selectively blocks so- 
dium current ,„ excitable cel., ft appears to block by occluding the extern* 
mou th f lh chilnnd porc> , by par . aI cn . nto >; ^J.^^t 

noncovalen ly to a turable carboxyl oxygen with a pK. of 5-6^ which b prob 
a. ly part ol the channel's selectivity filter. The selectivity filter is about 20-30% 
•n o the membrane electric field, and this causes TTX binding to be somewhat 
voyage dependent/' S TX competes with TTX, so th^ wo tL^I 

pelvl Th ,0 : S m tH . C ° UtSidC S0,Uti0n ' *" ^-tion appears to be com- 
pctmve '. The toxms have minimal effects on the channel ^ ting process as 

chanL k u l0Ck " USC - dc P cnd -t, ' as if their binding sL are 

thcTanneT * 1,10 C ° nf °"" ™d with'gating of 

Typieal nerve and skeletal muscle Na channels are blocked by TTX with 
1.1 stoichxometry w lt h K, of about 10 nM. However, cardiac Na Ihl^JZ 
^SSUB ^ with K, of a few „ M m^Lt ^ ^^^ Z Z 

m eh be a H-iT c ^ WM ^ *"* duc that thc cardiac Na ^annel 
r^ fs and t H ^ ^ haVe sh ™ that *e association 

low affinity of the carcbac channels is not unique. For example, Na vLnefa in 
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Figure I. STX dose-response curves for block of the cardiac (RHl) and the nerve 
(Brai»2a) Na channel CC-subuniis expressed in Xenopus oocytes. The STX dose-response 
curve for block of the cardiac isoform with Cys37+ (cardiac isoform numbering) mu- 
tated to a Plic is also shown. The currents were generated by trains of voltage steps 
using the iwo-pipctte voltage clamp method and the steady-state peak Na current is 
reported as a traction of the current in the absence of drug. 



the dorsal root ganglion belong to several populations, some with typical high 
TTX/STX affinity, some with low affinity, and some that are totally resistant. 121 



Sodium Channel Structure 

Na channels are scarce membrane proteins that are difficult to purify, so 
TTX and STX were very important markers for the Na channel during the long 
process of purification and amino acid analysis that was necessary to prepare for 
cloning of the protein. The first Na channel to be cloned was from the eel 
clcctroplax, and the Ot-subunit is a 2000 amino-acid protein. ,:,) Subsequendy, 
three Na channel a-subunits have been cloned from brain, two from skeletal 
muscle, and one from heart, all having 70-90% homology. ,4) The embryonic 
isoform in skeletal muscle, also manifested after denervation of adult muscle, 
appears to be identical to the cardiac isoform. These mammalian isoforms of the 
a-subunits have been successfully expressed in heterologous cells. The cardiac 
isoform has been shown to have the permeation, gating, and pharmacological 
properties of native cardiac Na channels. 1 5> A number of non-mammalian Na 
channel isoforms have been cloned, and other mammalian nerve and muscle 
isoforms have now been found. 14 '' Brain Na channels have two small J3 subunits, 
and skeletal muscle has one. 161 The presence or role of /3 subunits for the cardiac 
channel or the others has not yet been resolved. All isoforms of the a-subunits 
show a similar topology, determined by hydropathy analysis and the sidedness of 
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Figure 2. Likely topology of I he cardiac Na channel a subunii. Number of charges 
in the putative membrane spanning regions arc denoted. PKG site is denoted by 
Putative PKA nites (•) were taken as R/K X S/T or R X X S/T or R X S/T where 
X denotes any amino acid. Putative N-linkcd glycosylation sites (■) were defined by N 
X S/T where X is any amino acid except proline. (l ? mnv " with permission oi* the 
American Physiological Society) 



certain antibody, glycosylating and phosphorylation sites (Figure 2). There arc 
24 transmembrane segments, organized into four homologous sets of six seg- 
ments that appear to form domains. The fourth transmembrane segment of each 
domain has negatively charged amino acids in every third position, which led 
Numa and colleagues' : » to suggest that these segments are the channel's voltage 
sensor. 

The presence of four domains also led to the idea that the transmembrane 
path for ion transport was formed as a central pore that was surrounded by the 
four domains. Intuitive modeling of the molecule by Guy and colleagues 17 * sug- 
gested that the pore was partially lined by an infolding or the extracellular seg- 
ment between S5 and Sf> of each domain that they called SSI and SS2. In a 
brilliant set of mutational studies, Stuhmcr and colleagues , " ) proved that these 
segments were indeed part of the permeation path, and were also part of the 
binding site for TTX and STX. Two carboxyls in each of domains I and II, a 
carboxyl in domain IV, and a lysine in domain III were all critical for binding of 
TTX and STX, confirming Hillc\s proposal that the toxins bound to carboxyl 
oxygens in the pore mouth. Further, their mutation to neutral amino acids re- 
duced substantially the single channel conductance, as expected if the toxin site is 
in the permeation path. 

Isoform Differences In Guanidinium Toxin Binding 

Localization of the pore and toxin binding site, however, fa iled to explain 
thebyyi_ca£ diac affinity for the guanidinium toxins in comparison with hrain and 
^skeletal Na rhannrls. All of the residues identified as critical for toxin binding by 
the experiments of Terlau et al l8) are identical in the cardiac isoform. Only two 
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Figure 3. Depiction of the first domain or the Na channel a-subunit. The -seven 
amino acid SS2 sequence is shown and the differences between the cardiac isoform 
(RHI) and the brain and skeletal muscle isoform* (/il-Brainl!) isoforms arc indicated. 
Abbreviations for the amino acid residues are: C, Cys; D, Asp; £, Glu; F, Phc; N, Asn; 
Q, Gin; R, Arg; T, Thr, W. Trp; and Y ( Tyr. (From, 101 with permission or the Ameri- 
can AjBodiiiion lor the Advancement of Science.) 

residues differ between the isoforms in the four SS1-SS2 segments, a cysteine in 
place of a phenylalanine and an argininc in place or an asparagine (Figure 3). 
The cardiac channel also differs functionally in being especially sensitive to block 
by Cd J+ , and Cd i+ competes with STX binding/* Because cysteine is a logical 
amino acid for Cd* 2+ binding, Schild and Moczydlowski m) further suggested that 
this is the residue responsible for the difference in toxin binding. Within a short 
time, overwhelming evidence was found to prove that the residue in that location 
was critical: Satin ct al ,l,) replaced Cys in the cardiac isoform with Tyr or Phc and 
found an increase in toxin sensitivity to the levels seen in brain, Backx et aP> 
replaced the Tyr in the skeletal muscle isoform with Cys and reduced toxin 
sensitivity dramatically, and Hcincmann ct al 2l > replaced the Phc in the brain II 
isoform with Cys and reduced toxin sensitivity to the levels seen in heart. Change 
of the Arg/ Asn residue had only small effects on the toxin block, indicating that 
it had a minor role in the isoform difference. 



A Model of the Guanidinium Toxin Site 
At this point in our study of the toxin site the following information was 
available: 

1. Crystallographically determined structures of the rigid TTX and STX 
molecules. 

2. Knowledge of the toxin reactive sites from studies with toxin analogs. 

3. Identification of 6-8 residues in the channel pore of brain, skeletal and 
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cardiac isoforms that are required for toxin binding. 
4. Awareness of the specific residue responsible for isoform differences in 
binding. 

It was reasonable to ask if we could use this information set to infer the 
structure of the toxin binding site, using the toxins as structural templates of the 
sites. Presumably the toxins correspond to a key that fits securely into their locks 
by close association with channel residues, in order to achieve their high affinity 
and selectivity. If this effort to visualize the toxin binding site was successful, then 
it is possible that removal of the toxin from the site would reveal the general 
structure of the mouth of the channel pore. 

Our first task was to identify a possible secondary structure of the SSI-SS2 
regions such that the cridcal residues for toxin interaction in the four domains 
that have been identified by mutational studies could face a potential binding 
pocket. The motif for binding in domain I involved two adjacent residues-Asp 
and Cys/Tyr/Phc- and after an intervening Trp another critical residuc-Glu. 
This arrangement was unsuitable for OC-helical or ^-structures, but could occur in 
a reverse turn. We examined the probability of a /}-turn in the 5SI-SS2 segments 

I II 

R N T R 

L M v 

M W ^ ■ * I 1» 

383 Q G 941 



TTX/STX 



I L IV 
T G * _ D A 



T ^-M T 
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1714 B^tf 1422 

v v . 

IV III 

Figure 4. Structural motif for the TTX/STX binding site of the Na channel - an 
amiparallel arrangement of ^-hairpins of domains MV. Residues that interact with the 
toxins are marked by the arrows. The residue numbers are according ro the Brain2 
isoform, and are for the ones located in positions i + I of /3-turns. Relative orientation 
of /J-turns of the four hairpins is shown by flat arrows connecting the beta I and beta2 
strands, and it corresponds to an overhead view. {From 1 * with permission of the Bio- 
physical Society). 
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Figure 5. The IJpki.uW'o/./ar.l (rmKl.n.»»n l.i.ulmt? .silo model. TIX micrarw 
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and found a site that fit the Chou-Fasman rules Tor such a turn in each of the 
SS 1 -SS2 segments. Using software for molecular modeling, wc constructed lrom 
each segment /3-hairpins that were energetically stable** (Figure 4). After con- 
structing the toxin structures in the computer, wc could dock TTX onto the 
domain I /3-hairpin into an energetically satisfactory position and then dock the 
domain II 0-hairpin onto the complex. The resulting binding pattern (Figure 5) 
was a set of three salt bridges between the two carboxyb of domain I and the 
inner one of domain II to the guanidinium group of TTX. The two hydroxyl 
groups on C9 and CIO interacted by hydrogen bonds with the outer carboxyl of 
domain II. Each channel residue-toxin interaction was about - 10 kcal/mol un- 
der the conditions of the model construction. Superimposing the 7 8 
guanidinium of STX on the 1,2,3 guanidinium of TTX positioned the C12 
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hydroxyls of STX nearly exactly where the G9 and CIO hydroxyls of TTX were 
found. We therefore found that the same spatial arrangement of the domains I 
and II hairpins fit the STX molecule, but STX has a second guanidinium. This 
allowed docking of the domain IV hairpin onto the STX complex in a similar 
pattern, so that the domain IV carboxyl interacted with the second guanidinium 
of STX. The domain III hairpin was positioned symmetrically to allow modest 
non-bonding interacdons, forming a four-hairpin structure. This model fits the 
data of Terlau et al ,8) for the toxin binding site. 

If this four-hairpin model has any resemblance to the TTX/STX binding 
pocket, then the isoform differences seen as a consequence of the Cys/Tyr/Phe 
residues should be predicted by the model. Replacement ofCys with Tyr or Phe 
in the model positioned the aromatic ring to have a hydrophobic interaction with 
the toxin (Figure 5). The change in binding energy estimated by the modeling 
program was about - 5 kcal/mol, similar to the - 4 kcal/mol estimated experi- 
mentally by Satin ct al. m) 



Implications For Structure of the Channel 
Vestibule and Selectivity Filter 

Removal of the toxin from the binding site model revealed a l2A-dcep 
funnel-like structure about 12A on us side at the wide mouth, and narrowing to 
a 3 x 5A opening at the narrow mouth. This narrow opening corresponded to 
the size predicted by Hillc 21 * for the Na channel selectivity filter from the sixes of 
permeating organic cations. The opening was surrounded by four aligned resi- 
dues, one from each hairpin-Asp, Glu, Lys, and Ala. Heincmann ct al* IJ noticed 
the remarkable similarity of the SS I -SS2 regions of the Na and Ca channels, and 
similar alignment these four regions showed that the putative selectivity ring of 
the Na channel corresponded to four carboxyls in the Ca channel. They there- 
fore substituted carboxyls in these positions in domains III and IV of the brain II 
Na channel in place of the naturally occurring Lys and Ala, and found that they 
had created a mutant channel that was permeable to Ca- + . In their experiments 
the Lys to Glu mutation in domain III was the most effective. This dramatic 
demonstration of change in selectivity supported the idea that this ring of four 
amino acids in the Na channel constitutes part of its selectivity filter. 

We explored the idea further by allowing a model of the Na* ion hydrated 
with six water molecules to interact with the putative selectivity ring in our 
model. It readily interacted with the two inner carboxyls of domains I and II with 
an energy of about 100 kcal/mol-equivalent to the hydration energy of Na*. This 
interaction effectively dehydrated the ion, a step usually considered to be neces- 
sary for permeation. Ca 2 * interacted with the two carboxyls, but with insufficient 
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energy to dehydrate the ion. However, if the Lys of domain II was replaced with 
a carboxyl, the three carboxyls were positioned correctly to dehydrate Ca 2+ , as 
needed for its permeation, agreeing with the experiments of Heinemann et al. 24) 
K* ions interacted with the domain I and II carboxyls with the same energy as for 
Na*, so this analysis did not explain the 10:1 selectivity of the channel for Na* 
over K*. The model exhibits several characteristics necessary for the Na channel 
selectivity filter, but the model calculations cannot provide proof of its accuracy. 
More likely, the entire vestibule participates in the energetic steps that allow 
permeation. 

Several approaches to test and refine the Na channel vestibule-selectivity 
filter model arc being pursued. First, other Iigands that bind to the pore can be 
studied in mutated channels, including TTX/STX analogs and the /l-conotoxin 
peptides derived from snails. Second, location of component amino acid side 
chains can be explored by substitution with cysteines, with which divalent ions 
and organic sulfhydryl reagents may interact. Third, the nature and location of 
other residues involved in the permeation and/ or selectivity process can be found 
by detailed biophysical analysis of point mutations substituting residues of differ- 
ent, srze, hydrophobicity, and charge. Several of these efforts are now in 
progress.- 5 , -' 7) Similar efforts to model the K channel pore from toxin interaction 
and mutation of critical residues might be revealing, 281 because that channel has 
a much greater selectivity than the Na channel (about 100:1 for K + over Na + ). 
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